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Abstract 

We review the description of high energy diffraction from both the s and t channel 
viewpoints, and demonstrate their consistency. We emphasize the role played by s chan- 
nel unitarity and multi-Pomeron exchanges. We explain how these effects suppress hard 
diffractive processes. As examples, we describe the calculation of diffractive dijet pro- 
duction at the Tevatron and predict some of the rich diffractive phenomena accessible in 
proton-nuclear collisions at RHIC. 



1 Introduction 

A broad class of processes at high energies has properties analogous to the classical pattern of 
the diffraction of light. These are usually called diffractive processes. The classic example is the 
elastic scattering of hadrons on nuclei, which has an angular distribution with a series of minima 
and maxima, analogous to the diffraction of light on a black disc. In diffractive processes 

'To appear in the special issue of Acta Physica Polonica to celebrate the 65th Birthday of Professor Jan 
Kwiecihski. 
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the wave nature of particles is clearly revealed and the necessity of a quantum mechanical 
description in terms of amplitudes, and not probabilities, is evident. 

Unitarity relates the imaginary parts of forward elastic scattering amplitudes to total cross 
sections. Hence studies of elastic scattering and a to t for different targets and projectiles was 
traditionally one of important source of information on the strengths of interactions and their 
radii. In particular, data on the high-energy interactions of hadrons indicate that the radius 
of the strong interaction increases with energy, and indicate that the total cross sections will 
continue to increase in the energy region not explored at present accelerators. 

Processes of diffractive dissociation of colliding particles [1] provide new possibilities for 
the investigation of the dynamics of high-energy hadronic interactions. These processes are 
dominated by soft, nonpertubative aspects of QCD, and hence provide important information 
on the dynamics. Indeed, they provide a rich testing ground for models of soft interactions. 
In this paper we demonstrate that the cross sections of inelastic diffractive processes are very 
sensitive to nonlinear, unitarity effects. 

A new class of diffractive phenomena — hard diffraction — has been extensively studied both 
experimentally and theoretically in recent years. These processes include the total cross section 
and diffractive dissociation of a virtual photon at high energy at HERA, and the diffractive 
production of jets, W, Z bosons, heavy quarks, etc. in hadronic collisions. The latter processes 
become possible at very high (cm.) energies because large mass M diffractive states may 
be produced with M 2 <C s. Investigation of these processes gives new possibilities for the study 
of the interplay of soft and hard dynamics in QCD. An important question — the breaking of 
QCD factorization for hard reactions in diffractive processes — will be discussed below. 

Nowadays diffractive hard processes are attracting attention as a way of extending the 
physics programme at proton colliders, including novel ways of searching for New Physics, see, 
for example, [2] -[6] and references therein. An especially interesting process is the exclusive 
double diffractive production of a Higgs boson at the LHC; pp — > p + H + p, where the + signs 
denote the presence of large rapidity gaps. Clearly a careful treatment of both the soft and 
hard QCD effects is crucial for the reliability of the theoretical predictions for these diffractive 
processes. 

The outline of the article is as follows. In Section 2 we review diffraction phenomena from 
both the s and t channel viewpoints. In Section 2.3 we link these approaches together. We 
emphasize the convenience of working in terms of diffractive eigenstates, as originally proposed 
by Good and Walker [7]. In Section 3 we discuss how these eigenstates may be identified with 
the different parton configurations of the proton. Scattering through each of the eigenstates is 
suppressed (by unitarity or multi-Pomeron effects) by different amounts. In Section 4 we study 
some examples of hard diffractive processes. First we summarize a calculation of diffractive 
dijet production at the Tevatron, which shows that the estimate of the overall suppression of 
the cross section is in agreement with the data. In Section 4.2 we turn to diffractive proton- 
nuclear processes. We present predictions for a range of hard diffractive processes which are 
accessible at RHIC. Finally, Section 5 contains some brief conclusions. 
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Since this article is an attempt at a short review 1 of high energy diffractive phenomena, we 
have only provided a few references to help the reader. Hence our reference list is incomplete 
and we apologize to those authors whose work has not been cited. Section 4 contains new 
results for diffractive proton-nuclear processes. 



2 Theoretical approaches to diffraction at high energies 

The investigation of the diffractive scattering of hadrons gives important information on the 
structure of hadrons and on their interaction mechanisms. Diffractive processes may be studied 
from either an s channel or t channel viewpoint. We will discuss these approaches in turn, and 
then show that they are complementary to, and consistent with, each other. 



2.1 Diffraction from an s channel viewpoint 



Unitarity plays a pivotal role in diffractive processes. The total cross section is intimately 
related to the elastic scattering amplitude and the scattering into inelastic final states via 
s channel unitarity, SS^ = I, or 



disc T = T — T t = iT^T 



(1) 



with S — I + iT. If we were to focus, for example, on elastic unitarity, then disc would simply 
denote the discontinuity of T across the two-particle s channel cut. At high energies, s channel 
unitarity relation is diagonal in the impact parameter, b, basis, such that 



2ImT cl (s, b) = \T el (s, b)\ 2 + G M (s, b) 



with 

Ctot 
Cinel 

The general solution of (2) is 



2 J d 2 blmT cl (s,b) 
J d 2 b\T cl (s,b)\ 2 

2lmT cl (s,b)-\T cl (s,b)\ 2 



d 2 b 



Tel 
G^inel 



i(l - r]e 2iS ) 

i-v 2 , 



(2) 

(3) 
(4) 
(5) 



(6) 
(7) 



where i](s,b) 2 = exp(— Q(s, b)) is the probability that no inelastic scattering occurs. Q > is 
called the opacity (optical density) or eikonal 2 . 

x Other reviews of diffraction can be found, for example, in [8]-[14]. 
2 Sometimes 0/2 is called the eikonal. 
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The well known example of scattering by a black disc, with r] = for b < R, gives a e \ = 
Cinei — 71 R 2 and o"tot = 27ri? 2 . In general, we see that the absorption of the initial wave due 
to the existence of many inelastic channels leads, via s channel unitarity, to diffractive elastic 
scattering. 

So much for elastic diffraction. Now we turn to inelastic diffraction, which is a consequence 
of the internal structure of hadrons. This is simplest to describe at high energies, where the 
lifetimes of the hadronic fluctuations are large, r ~ E/m 2 , and during these time intervals the 
corresponding Fock states can be considered as 'frozen'. Each hadronic constituent can undergo 
scattering and thus destroy the coherence of the fluctuations. As a consequence, the outgoing 
superposition of states will be different from the incident particle, and will most likely contain 
multiparticle states, so we will have inelastic, as well as elastic, diffraction. 

To discuss inelastic diffraction, it is convenient to follow Good and Walker [7], and to 
introduce states <p k which diagonalize the T matrix. Such eigenstates only undergo elastic 
scattering. Since there are no off-diagonal transitions 

<^|T|0 fc > = O iotj^k (8) 

a state k cannot diffractively dissociate in a state j. We have noted that this is not true for 
hadronic states due to their internal structure. One way of proceeding is to enlarge the set 
of intermediate states, from just the single elastic channel, and to introduce a multichannel 
eikonal. We will consider such an example below, but first let us express the cross section 
in terms of the probabilities Fk of the hadronic process proceeding via the various diffractive 
eigenstates <p k . 

Let us denote the orthogonal matrix which diagonalizes ImT by a, so that 

ImT = aFa T with {<j>j\F\<j> k ) = F k 5 jk . (9) 
Now consider the diffractive dissociation of an arbitrary incoming state 

K) = S aik l^ fc )- ( 10 ) 

k 

The elastic scattering amplitude for this state satisfies 

(i\lmT\i) = ]T \a lk \ 2 F k = (F), (11) 
k 

where F k = {<p k \F\<j) k ) and where the brackets of (F) mean that we take the average of F over 
the initial probability distribution of diffractive eigenstates. After the diffractive scattering 
described by T^, the final state |/) will, in general, be a different superposition of eigenstates 
from that of \i), which was shown in (10). At high energies we may neglect the real parts of 
the diffractive amplitudes. Then, for cross sections at a given impact parameter b, we have 

^ = 2Im(i|r|t> = 2 £ \a tk \ 2 F k = 2(F) 
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i\T\i)\ z = V \a tk \ 2 F k = (F) (12) 



d 2 b 



It follows that the cross section for the single diffractive dissociation of a proton, 

^ = (F>) - (F}\ (13) 

is given by the statistical dispersion in the absorption probabilities of the diffractive eigenstates. 
Here the average is taken over the components k of the incoming proton which dissociates. If 
the averages are taken over the components of both of the incoming particles, then (13) is the 
sum of the cross section for single and double dissociation. 

Note that if all the components <pk of the incoming diffractive state \i) were absorbed equally 
then the diffracted superposition would be proportional to the incident one and the inelastic 
diffraction would be zero. Thus if, at very high energies, the amplitudes at small impact 
parameters are equal to the black disk limit, F k — 1, then diffractive production will be equal 
to zero in this impact parameter domain and so will only occur in the peripheral b region. 
Such behaviour already takes place in pp (and pp) interactions at Tevatron energies. Hence 
the impact parameter structure of inelastic and elastic diffraction is drastically different in the 
presence of strong s channel unitarity effects. 

Under the assumption that amplitudes F k at high energies can not exceed the black disk 
limit, Fk < 1, equations (12) lead to the following bound 

^o"c1 + sd 1 da^t /,.\ 
d?b -2^Pb' [ ' 

known as the Pumplin bound [15]. 

A simple realistic application of the above framework is the interaction of a highly virtual 
photon with a nucleus. For recent discussions and references see, for example, [16, 17]. The 
wave function for the photon to fluctuate into a qq pair may be written ip(r, a), where f is the 
separation of the q and q in the transverse plane and a, 1 — a are the longitudinal momentum 
fractions of the photon momentum carried by the q and q. For small r the cross section for the 
qq pair to interact with a nucleon behaves as a ~ a^r 2 . Thus, at fixed impact parameter b, the 
opacity 

Q = aT A (b) oc r 2 (15) 

where 

roo 

T A (b) = / dzp(b,z) (16) 

and p(b, z) is the density of nucleons in the nucleus (/ T A (b) d 2 b = A). T A {b) is often called the 
nuclear density per unit area or the optical thickness of the nucleus. 
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At high energy, due to Lorentz time dilation, the lifetime of each component with fixed r 
is much longer than the interaction time of the qq pair inside the nucleus A. Hence the 7* 'A 
amplitude is of the form 3 

(1 - e- n ' 2 ) = J d 2 rda*fj*{r,a) [1 - exp{-a{r)T A {b)/2)] i/>(r,a). (17) 

On the other hand, at low energy, where the lifetime of the qq fluctuations is small, we have 
to average a(r) already in the exponent, since each nucleon interacts with a different 7* — > qq 
fluctuation. Therefore the low energy amplitude is 

l-exp(-(fi)/2) = l-exp(-( ( 7)r A /2). (18) 

To simplify the discussion we, here, neglect the real part of the elastic amplitude; that is we 
assume in (6) that 5 <C 1. Modulo this assumption, (18) is exactly the Glauber formula, which 
was originally derived for an incoming particle of energy E < 1 GeV, which therefore has 
enough time to return to its equilibrium state between interactions within the nucleus. It is 
this normal mix of Fock states of the incident particle which leads to taking (a) in the exponent. 

Returning to the high energy case, we see that as the separation r is frozen, each value 
of r corresponds to a diffractive eigenstate (in the notation of (10)) with cross section 
Ck = &(f)- If there is a large probability of inelastic scattering, that is crT > 1, then the 
difference between (17) and (18) is significant, and leads to a large diffraction dissociation cross 
section usd of (13). 

In the above example we are dealing with a continuous set of diffractive eigenstates \<pk), 
each characterised by a different value of r. Another example is to consider just two diffractive 
channels [20]-[23] (say, p, N*), and assume, for simplicity, that the elastic scattering amplitudes 
for these two channels are equal. Then the T matrix has the form 

Im T = 1 - e' n '\ (19) 

where the eikonal matrix Q has elements 

n f j v = Q uj fi uj ri ,. (20) 

The individual u> matrices, which correspond to transitions from the two incoming hadrons, 
each have the form 

«=(;?)■ « 2i) 

The parameter 7(5, b) determines the ratio of the inelastic to elastic transitions. The overall 
coupling fl Q is also a function of the energy y/s and the impact parameter b. It is assumed here 

3 This result was first proposed within QCD in Refs. [18] and [19], but, long before, an analogous expression 
was originally given [1] for deuteron scattering with proton and neutron constituents in the place of the qq 
fluctuation of the photon. 
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that the diagonal elements for both channels are equal. This leads to a simplification of the 
formulas. 

With the above form of u, the diffractive eigenstates are 

= ^ (b) + \N*)) , I0 2 ) = ^= (b) - . (22) 
In this basis, the eikonal has the diagonal form 

^m"n' = ^0 d mn d m > n >, (23) 

where m, n = 4>ii4>2 an d 

In the case where 7 is close to unity, 7 = 1 — e, one of the eigenvalues is small. As was mentioned 
above, in QCD the diagonal states are related to the colourless qq states with definite transverse 
size, r. Small eigenvalues correspond to the states of small size (a ~ r 2 ). Thus the s channel 
view of diffraction is convenient to incorporate s channel unitarity. However it needs extra 
dynamical input to predict the s and b dependences of the different diffractive processes. 



2.2 Diffraction from a t channel viewpoint 

The t channel approach is based on the Regge model for diffractive processes. In this approach 
diffractive processes are mediated by the exchange of a Pomeron (P) - the leading Regge pole 
with vacuum quantum numbers (Fig. 1). The Pomeron plays the role of an exchanged 'particle', 
and gives factorizable contributions to scattering amplitudes. 

Diffractive processes (Fig. 1) are characterized by a large rapidity gap between groups of 
produced particles. For example, for the (single) diffractive dissociation of particle 2 into a 
system of mass M, the rapidity gap between particle 1' and the remaining hadrons is 

^ = = (25) 

where xp = 1 — x and x = Pvl/ (pi'L)max is the Feynman x variable of particle 1'. In general, 
the masses, Mj, of the diffractively excited states, produced in high energy ^/s collisions, can 
be large. The only condition for diffractive dissociation is Mf s. 

In the Regge pole model, the cross section for the inclusive single diffractive dissociation 
process of Fig. 1(a) can be written in the form 

= ^l«*(AM)l^(M>, t ) (26) 

where t is the (square) of the 4-momentum transfer. The Green's function (propagator) of the 
Pomeron is 

G P (Ay,t) = Sexp[(a P (t) - l)Ay], (27) 
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Figure 1: Single and double diffractive dissociation with a large rapidity gap represented by 
Pomeron exchange. 



where 

s = 1 + exp{-inap(t)) 
sin 7r ctp(t) 

is the signature factor. The quantity a^^M 2 ,t) can be considered as the Pomeron-particle 2 
total interaction cross section [24]. This quantity is not directly observable. It is defined by 
its relation to the diffraction production cross section, (26). This definition is useful, however, 
because at large M 2 , this cross section has the same Regge behavior as the usual total cross 
sections 

/ M 2\ «*(0)-l 

41(M 2 ,t) = £<&(0)r PP (*) — (29) 
k \ s o J 

where the r PP {t) is the triple- Reggeon vertex, which describes the coupling of two Pomerons 
to the Reggeon k. 

Thus, in this region where the dissociating system has a large mass M, satisfying s ^> 
M 2 3> m 2 , the inclusive diffractive cross section is described by the triple- Regge diagrams of 
Fig. 2. Therefore, on inserting (29) into (26), we have 

^£r iW (i-.p~(^) . (30) 

The Pomeron and / couplings to the proton, and the triple-Regge vertices r PP ,r PP , have 
been determined from analyses of experimental data on the diffractive production of particles 
in hadronic collisions (see the review in [8]). These couplings specify the factors Vj-it) in (30). 

In impact-parameter space, Regge amplitudes have a gaussian form at asymptotic energies. 
On the other hand, for inelastic diffraction we expect a peripheral form, which follows from 
unitarity in the s channel picture. Note also that, since the Pomeron has intercept ap(0) > 1, 
the cross sections of diffractive processes, (30), increase with energy faster than the total cross 
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Figure 2: Triple- Regge diagrams with k = P or / exchange describing high energy (M 2 3> m 2 ) 
Pomeron-particle 2 scattering. 

sections, and thus lead to a problem with unitarity. However, in Regge theory it is necessary 
to take into account, not only Regge poles, but also Regge cuts [25, 26], which correspond to 
the exchange of several Regge poles in the t channel. These contributions restore the unitarity 
of the theory. For inelastic diffraction they lead to a peripheral form of the impact-parameter 
distributions and yield the appropriate energy dependence of the corresponding cross sections. 

How is this realized in Regge theory? The explanation is based on a technique to evaluate 
Reggeon diagrams, introduced by Gribov [27], which allows the calculation of contributions of 
multi-Pomeron cuts in terms of Pomeron exchanges in the amplitudes of diffractive processes. 

2.3 Compatibility of the s and t channel viewpoints: 

Gribov's Reggeon calculus and the AGK cutting rules 

Gribov's technique uses unitarity and analyticity of Reggeon-particle amplitudes, which fol- 
lows from an analysis of Feynman diagrams [27]. For example the amplitude of two- Pomeron 
exchange in the t channel (Fig. 3a) can be written as a sum over all diffractive intermediate 
states in the s channel (Fig. 3b). This result can be extended to multi-Pomeron exchanges 
in the t channel. For the single (s) channel case, the summation of all elastic rescatterings 
leads to the well known eikonal formula with eikonal given by Fourier transform of the single- 
Pomeron-exchange diagram. For several diffractive channels the resulting amplitudes can be 
written in the matrix eikonal form of the type shown in (19) and (20). In this way a connection 
between the s channel view of diffraction and Regge theory is established. Note that only for 
weak coupling can the Pomeron cuts be neglected, and just the first term in the expansion be 
retained. Another powerful tool for the calculation of, not only diffractive, but inelastic pro- 
cesses, which contribute to total cross section, is the AGK (Abramovsky, Gribov, Kancheli [28]) 
cutting rules. By applying these rules, it is possible to show the selfconsistency of the approach, 
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2 2 2 2 2 2 

(a) (b) 

Figure 3: Two-Pomeron exchange in the t channel expressed as the sum over all diffractive 
intermediate states in the s channel. The crosses indicate that the particles are on the mass 
shell. 



which was absent in the pure Regge pole model. Consider a diagram where elastic scattering 
is mediated by the exchange of k Pomerons. The AGK cutting rules specify the coefficients 
c]J arising when n of these Pomerons are cut. Recall that these Pomeron cut discontinuities 
give the corresponding inelastic contributions to a to t- The terms with n = correspond to the 
diffractive cutting of the diagram (that is the cut is between the Pomeron exchanges, and not 
through the Pomerons themselves). These n = terms have the coefficients 

c n=0 = gfc-l _ L ( 31 ) 

For two-Pomeron exchange, k = 2, the coefficients are +1, —4, +2 according to whether n = 0, 
1 or 2 Pomerons are cut, respectively. 

We illustrate the compatibility between the s channel view of diffraction and the t channel 
approach with multi-Pomeron cuts, using a simple example with two diffractive channels. In 
this case the T matrix has the form (19)— (21) with function fl given by 

n °~ 4^B 6 - I^ eXP (, Aln UJ~45j' (32) 

where B is the slope of the Pomeron amplitude, 

B = ±B + a' P ln(s/s ), (33) 

g£p is the coupling of the Pomeron to the proton, s = 1 GeV 2 and A = Op(0) — 1. In practice, 
the value of 7 is about 0.4. Thus as a first approximation we may neglect higher powers of 7 
in (19). The eikonal function Q increases with energy as (s/so) A , and at very high energies it 
follows, from (32) and (33), that Q ^> 1 for values of the impact parameter satisfying 4 

b 2 <R 2 = Aa'pA In 2 (s/s ) . (34) 
4 We have neglected the slowly varying term a' P ln(s/so) m [(9p P ) 2 / ^ na 'p hi(s/so)] in R 2 . 
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As a consequence, the elastic and inelastic cross sections have strikingly different forms in impact 
parameter space. The elastic cross section is limited by the black disc bound, a e \(b,s) = 1, 
up to b ~ R, after which it decreases rapidly. On the other hand, the inelastic cross section 
depends on the off-diagonal elements in (20) and (21), 

a 12 (b,s) ~( 7 fio) 2 e- n °, (35) 

and so is only non-zero in the peripheral region of the interaction, b ~ R. This example 
illustrates the general difference in the suppression of the elastic and inelastic channels due to 
multi-Pomeron exchanges, in agreement with s channel unitarity. 

Let us illustrate the self consistency of the Gribov technique for the evaluation of multi- 
Pomeron contributions, and the AGK cutting rules, in the simplest single-channel example. In 
this case the imaginary part of elastic amplitude in b space has the eikonal form 



I m T cl (s,b) = J2 { - 1 ^{Y) =1 - ex P(-^o/2). (36) 



By applying the AGK cutting rules [28] to each term of the series, it is possible to obtain for 
the diffractive (elastic) cross section [29] 

a A (s,b) = E ^p (y)* 2 ^*" 1 " 1) = (1 " exp(-fio/2)) 2 (37) 
and for inelastic cross section 

<w(a, 6) = E = (1 - exp(-no)), (38) 



n=l n - 

where n is the number of cut Pomerons (with any number of uncut ones). It has a simple 
probabilistic interpretation. The second expression for cr inc i represents the whole probability, 1, 
minus the probability exp(— Qq) to have no inelastic interaction, whereas in the first expression 
each term Qq/uI represents the probability of n inelastic interactions (where n\ accounts for the 
identity of the interactions) multiplied by exp(— Q ) which guarantees that there are no further 
inelastic interactions. 

/,From (36)-(38) we see that, in the single channel case, we have 

a el (s,b) = (ImT el ( S ,fo)) 2 (39) 

a c i(s, b) + a- mc i(s, b) = a tot (s, b) = 21m T cl (s, 6), (40) 

as indeed it must be. Note that this self consistency could not be achieved if the series in 
multi-Pomeron rescatterings would be cut at some finite term. The consistency with s channel 
unitarity can be proven for an arbitrary multichannel case. 
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It is also instructive to consider the situation with two Pomeron poles: soft P s and hard P^. 
In the single channel case, a summation of s channel rescatterings leads again to the eikonal 
formula [30] with 

n( s , b) = n s ( s , b) + n h ( s , b) (41) 

where the Qi(s, b) are the Fourier transforms of the i-Pomeron contributions. The hard Pomeron 
is usually related to the cross section of minijet production [31, 30], which rapidly increases 
with energy. It is interesting to note that the unitarization via eikonal formula described above 
leads to a rather weak influence of the hard component on the behaviour of the total cross 
sections, because this component becomes important only at energies when the soft component 
is large and the amplitudes at small impact parameters are close to the black disk limit. From 
the AGK cutting rules it is easy to prove that in this model, the cross section of a hard process 
(with any number k > 1 of hard interactions plus any number of soft inelastic interactions) is 
described by the eikonal expression with f2 = fi h . This is the so-called self absorption theorem 
for processes with a given criterion [32]. In this case the criterion is the existence of at least 
one minijet. ^From the above expression it is clear that such processes will dominate at very 
high energies, and that the mean number of minijets will be determined by Q^. 



3 Diffractive eigenstates and parton configurations 

The partonic picture of hadronic fluctuations mentioned above gives a new insight to pro- 
cesses of diffraction dissociation of hadrons. This picture naturally appears in QCD, where 
each hadron can be represented as a Fock state vector in terms of the quark and gluon de- 
grees of freedom. At very high energies, these fluctuations have large lifetimes and have small 
variations during the time of the interaction. Thus configurations made of definite numbers of 
quarks, antiquarks and gluons with definite transverse coordinates can be considered as natural 
candidates for the eigenstates \<pk) of diffraction. Partonic models of diffraction were originally 
introduced in refs. [33, 34] and are now widely used within a QCD framework (for examples and 
references see [35]). In QCD, colourless configurations of small transverse size have small total 
interaction cross sections. Existence of such configurations leads to the phenomenon known 
as 'color transparency' [18, 19], which, in turn, leads to interesting effects in the diffractive 
interactions of hadrons and photons with nuclei. Existence of the diffractive eigenstates with 
substantially different interaction cross sections leads, according to (13), to large total cross 
sections of diffractive processes at high energies, which are close to the Pumplin bound (14). 
The distribution, P(<7k), of cross sections a k of the diffractive eigenstates <pk, has been deter- 
mined from experimental data on diffractive processes for pions and protons in Refs. [36]. In 
the simplest 2-channel model, discussed in Section 2.1, the broad distribution P(ak) is ap- 
proximated by 2 states with substantially different values of In ref. [37] these states were 
related to the distributions of quarks and gluons with substantially different values of Bjorken 
x (see also [38]). In Ref. [37] two models were considered. In the first model (model A) it was 
supposed that the valence quarks correspond to the small size component, while the sea quarks 
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and gluons make up the large size component. In the second model (model B) the weights of 
the two components varied with x in such a way that the small size component dominated at 
x ~ 1, while the large size component was dominant at small x. 

Another interesting aspect of the partonic picture is the existence of partonic fluctuations 
with large masses, M ^> m^, which correspond to multigluon configurations. In perturbative 
QCD, these are related to the BFKL Pomeron [39]. The elastic scattering of such states leads 
to the diffractive production of states of large mass, which in the Regge language corresponds 
to the triple-Pomeron interaction. The average mass of such states increases with energy, and 
at very high energies their role in inelastic diffractive processes becomes very important. 

It was emphasized above, that at very high energies when elastic scattering becomes close 
to the black disk limit, there is a strong influence of different diffractive eigenstates on the 
properties of diffractive processes in the physical (hadronic) basis. In particular, there is a 
strong reduction of the cross section of large mass diffraction due to unitarity (multi- Pomeron 
rescatterings) effects. Below we will consider these effects for the diffractive production of large 
mass states, which contain a hard subprocess: the so-called hard diffractive processes. 



4 Hard diffraction in pp (pp) and pA collisions 

The hard diffractive dissociation of hadrons or nuclei provides new possibilities for the inves- 
tigation of diffraction dynamics. We illustrate this by discussing some typical reactions below. 
First, we describe a study of diffractive dijet production at the Tevatron, and then we proceed 
to make predictions for typical hard diffractive processes in proton-nuclear collisions that are 
accessible at RHIC. The possible relevance to diffraction at the LHC is mentioned. 



4.1 Diffractive dijet production in pp collisions 

Consider first the situation for dijet production in pp (pp) collisions. It is important that the 
single Pomeron exchange diagram of Fig. 4(a) can be calculated using QCD factorization for 
the hard processes, together with the distributions of partons in the proton and the Pomeron, 
denoted by ff and ff respectively. The latter 'diffractive' distributions can be extracted from 
an analysis of hard diffraction in deep inelastic scattering (Fig. 4(b)). Thus the cross section 
for diffractive dijet production, Fig. 4(a), may be written 

= E/ F P (x P ,t)f l p (P,E 2 T )f k (x 1 ,E 2 T )a lk dpdx 1 (42) 

where is the cross section for dijet production by partons with longitudinal momentum 
fractions x\ and (3 of the proton and Pomeron respectively, and Et is the transverse energy of 
the jets. Fp is the Pomeron "flux factor", which, according to (26), can be written in the form 

Mx P ,t)= ^^rJl^ . (43) 



(lQit)x P 
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(a) (b) 

Figure 4: (a) The Born diagram for diffractive dijet production at the Tevatron, and (b) the con- 
tribution to diffractive production at HERA driven by the gluonic component of the Pomeron. 
Notice that the lower parts of the diagrams are the same. 

This approximation corresponds to the Ingelman-Schlein conjecture [40]. The existence of a 
hard scale provides the normalization of the Pomeron term. 

In the previous sections we have emphasized the importance of multi-Pomeron contribu- 
tions for diffractive particle production. For diffractive dijet production they correspond to the 
diagrams shown in Fig. 5. These diagrams lead to a strong violation of both Regge and hard 




Figure 5: Multi-Pomeron contributions to diffractive dijet production at the Tevatron. 

factorization. Experimental data of CDF Collaboration at the Tevatron [41] clearly indicate 
the breaking of hard factorization in dijet production — the measured dijet cross section is an 
order of magnitude smaller than that predicted from (42) using HERA data for ff. It is even 



14 



more striking that the (3 distribution for the produced dijets increases much faster as (3 — * 0, as 
compared with the expected behaviour of the partonic distributions. Both experimental obser- 
vations can be reproduced in the model which takes into account the suppression of diffractive 
production at high energies due to multi-Pomeron exchanges. An essential feature of the model, 
which allows the reconciliation of the observed (3 dependence with that expected from the par- 
tonic (mostly gluonic) distributions, is the relation between the size of the partonic fluctuation 
and the type (or the value of x) of the active parton. The suppression factor S p in this model 
can be written in the form 5 

J2 d 2 b\a pn \ 2 \M n \ 2 exp(-Q„(s, &)) 

S P = ~ —7 , (44) 

]T d 2 b\a pn \ 2 \M n \ 2 

n J 

where \a pn \ 2 is the probability of finding the partonic diffractive eigenstate n = \<p n ) in the 
proton \p), see (10), and |.M„| 2 is the probability of producing the dijet system from the 
eigenstate n. It can depend on the type of the active parton i and the value of its momentum, 

— * 

Xi and k±i. 

In Ref. [37] two eigenchannels were considered, n — 1,2. It was argued that the channel 
with the smaller cross section (denoted the S component) corresponds to valence quarks with 
x ~ 1, while the channel with the larger cross section is due to sea quarks and gluons and is 
concentrated at smaller values of x. This was denoted model A of the diffractive eigenstates. 
Of course, the model is oversimplified. There will be part of the valence component with large 
size, while on the other hand the gluons and sea quarks contribute to the small size component. 
An alternative model, model B, was introduced, in which the partonic distributions 

fi(x,E%.) = (Pf + P?)Mx,E%) (45) 

with i = valence, sea, glue; where the large and small size projection operators have the forms 

Pf = (1 - x )<h&\ P? = \-Pt (46) 

where the were chosen so that each component carries one half of the whole nucleon energy 
and one half of the number of valence quarks, see [37]. 

The functions Q n (s,b) in (44) have been parameterized in the form (see (23)) 

Q L = (1 + 7 K, n s = (1-7)^0, (47) 

with 7 = 0.4 [42]. In practice we use a more realistic parameterization of the Q n , determined 
from the global description of total, elastic and soft diffractive production data [42]. In addition 
to the two-channel eikonal allowing for low mass diffraction, this analysis incorporated pion-loop 

5 We use the subscript p to distinguish the suppression S p in pp (or pp) hard diffractive collisions from the 
extra suppression Sa in p-nuclear collisions to be discussed in Sections 4.2 and 4.3. 



15 



insertions in the Pomeron trajectory, and high mass single- and double-diffractive dissociation 
via s channel iterations of diagrams containing the triple-Pomeron interaction. 

Interestingly, the order of magnitude suppression and the (5 dependence of the prediction 
of the CDF dijet was well reproduced by both models A and B for the diffractive eigenstate, 
see Fig. 6, which was taken from Ref. [37]. Note that, when the new fits to HI diffractive 




Figure 6: The four lower curves are the predictions for diffractive dijet production at the 
Tevatron obtained from two alternative sets of 'HERA' diffractive parton distributions I and II 
using models A (continuous curves) and B (dashed curves) to calculate the suppression factor 
S p of (44). Note that without the suppression factor the upper two curves in the plot would be 
obtained. The Tevatron dijet data (shown by the data points and shaded band) demonstrate 
the importance of including the shadowing corrections. The figure is taken from [37]. 



data [43] are used in the approach of Ref. [37], even better agreement with CDF dijet data is 
achieved [44]. 

A further check of this general approach has been performed in Ref. [45]. There it was shown 
that the experimentally observed [46] breakdown of factorization in the ratio of the yields of dijet 
production in single diffractive and double-Pomeron-exchange processes is in good agreement 
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with the expectations of the suppression factors obtained from the above-mentioned global 
analysis of diffractive phenomena [42]. 



4.2 Diffraction in proton— nuclear collisions 

The suppression factors, analogous to (44), are much stronger for diffractive p-nuclear processes 
than for pp (or pp) collisions. They have a much richer structure. Their magnitudes, their A 
and kinematic dependences are much more sensitive to the absorption cross sections of the 
proton partonic configurations of different transverse size, see, for example, [47] . Moreover, for 
collisions involving nuclei, A, there are more possibilities of diffractive dissociation. We will 
consider two types of proton-nuclear hard diffractive processes: the incoherent and coherent 
production of a massive system h (accompanied by other particles, denoted by X) 

pA -> Xh + Y A , (48) 

pA -> Xh + A, (49) 
as shown in Fig. 7. The observable system h could be a dijet, a Drell-Yan pair or a heavy qq 




Figure 7: Multi-Pomeron exchange contributions to the incoherent and coherent diffractive 
production of a massive system h in high energy proton-nuclear (pA) collisions. In (a) the 
nucleus A has dissociated into a system Ya of nucleons, whereas in (b) it remains intact. The 
dots are to indicate any number of exchanged Pomerons, and/or interactions with any number 
of nucleons. 



pair. For incoherent production, (48), the nucleus dissociates into a system of nucleons labelled 
Ya- The momentum transfer is t ~ 1/Rp, where R p is the radius of the proton. On the other 
hand, in process (49) we have the coherent dissociation of the proton, while leaving the nucleus 
intact with very small momentum transfer, t ~ 1/R% where Ra is the nuclear radius. 
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In the Born approximation for the interaction with a single nucleon, the cross section, 
integrated over t, of the second process is 



~ £2/3 

t=o 



r a 2 

' ' = CM 4 / 3 , (50) 



while the cross section of the incoherent process behaves as a(pA — > Xh + Y A ) ~ A. These 
A-dependences are drastically modified by multiple rescattering (multi-Pomeron exchanges), 
shown in Fig. 7. 

The cross section for the incoherent hard diffraction on nuclei, process (48), in the Gribov- 
Glauber approach [48] is described by an analogous formula to that for pp collisions, 



a(pA — > Xh + Y A ) = £ / d 2 b \a pn \ 2 a n {pp - h)T A (b) exp(-< bs (s)T A (6)), (51) 

where, as before, the sum is over the partonic diffractive eigenstates of the incoming proton, and 
\a pn \ 2 is the probability of finding eigenstate n = \<p n ) in the proton \p), see (10); a n (pp — > h) 
is the corresponding cross section in proton-proton 6 collisions for scattering through the n th 
diffractive eigenstate, and 7a (6) is the nuclear profile function (16). Here a™ bs (s) = <r™ ot — 
a"^, where the cross sections are for proton-nucleon scattering at energy y/s through the n th 
eigenstate, since elastic scattering does not change the kinematic structure of the process. Thus 
the nuclear suppression factor (sometimes called nuclear transparency) is 

Q incoh _ °{V A -±Xh + Y A ) , ko , 
bA ~ Aa(pp-,Xh) ' (52j 

where the denominator a(pp — > Xh) is the corresponding cross section for the diffractive 
production of the hard (massive) system h in proton-nucleon collisions. Thus SA Coh is the extra 
suppression due solely to nuclear effects. 

Coherent production on nuclei, process (49), is possible when the mass of the diffractively 
produced system, M(hX), satisfies 

M 2 /s < l/m p R A , (53) 
where ^Js is the proton-nucleon cm. energy. Then the amplitude in impact parameter space 

M coh (pA ^Xh + A) = £ a pn M n (pp - h)T A (b) exp (-^T A (b)j , (54) 

where again the sum is over the diffractive eigenstates of the incoming proton. Here a^ bs = a™ ot , 
since even elastic scattering with t > 1/R A wm destroy the nucleus. Hence the total cross 
section of coherent diffractive production of h is 



is 



a coh (pA ^Xh + A) = J2 I d2b Kn\ 2 4n ^ 



da n (pp — > h) 



T 2 A (b)eM-<tT A (b)). (55) 

t=o 



a The Pomeron couples equally to the proton and neutron constituents of the nucleus. 
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We present below the results for the nuclear suppression factor for the coherent diffraction 
production of a system h in the form 



Again this is the extra suppression arising solely from nuclear effects. 

To be realistic, we present results for proton-nuclear (pA) collisions in the RHIC energy 
regime, yt~s ~ 300 GeV, where yHi refers to the corresponding proton-nucleon cm. energy. We 
use the same model that was employed to describe hard diffraction in pp collisions [42]. We 
show results where the system h is chosen to be Drell-Yan lepton pairs (which originate from 
qq annihilations), and also where h is chosen to be either cc or bb pairs (which originate from 
gg fusion). To be specific, we use a factorization scale /i 2 = Mp air where for Drell-Yan or cc 
production the mass of the pair is chosen to be M pair = 5 GeV, while for beauty production 
we take M b i = 13 GeV. The nuclear profile function, Xa(6) of (16), was calculated using the 
standard Woods-Saxon form for p(r) with the parameters determined from the electromagnetic 
form factors of nuclei [49]. In Figs. 8, 9 we show the nuclear suppression factors (52) and (56) 
as a function of the fraction of the longitudinal momentum of the incoming proton carried by 
the system h (= , cc or bb). Figs. 8 and 9 correspond to nuclei with A = 64 and A = 197 
respectively. We see that, in general, the suppression due to nuclear effects is approximately 
a factor of 10 for pA — > Xh + Y A and a factor of 100 for pA — > Xh + A. The dashed lines in 
Figs. 8 and 9 correspond to the predictions of the one-channel eikonal model, which obviously 
cannot represent the different partonic configurations participating in the process. It therefore 
does not depend on the specific diffractively produced system h or the momentum fraction x 
that it carries. The other curves are the predictions of the two different two-channel eikonals 
introduced in Section 4.1 — the simple model A, which predicts the dotted curves, in which the 
valence quarks and the gluon + sea quarks are respectively associated with the small and large 
size diffractive eigenstates, |02) and \<pi) of (47); and model B, which predicts the continuous 
curves, with a more realistic partonic composition of the diffractive eigenstates. 

For models A and B, we see that, in general, the nuclear suppression depends on x. For 
large x, where the small-size, low a a b s component dominates, the suppression is less than that 
expected in the pure Glauber one-channel case with exp(— {(j)Ta)- On the other hand, for 
smaller x, where more contribution comes from the large-size, higher (X a b s component, there 
is much more suppression. However, for charm (and beauty) production, model A predicts 
that the suppression is independent of x. The reason is that the gluons have all been assigned 
to a single component; moreover, as this component has a a b s > (c), the suppression is larger 
than the single-channel prediction. In the more realistic model, model B, where the gluon is 
distributed between both eigenstates, the suppression is less and the pA cross section larger, 
and moreover depends oris — the large x gluons are concentrated in the diffractive eigenstate 
with the small <r a b s - On the other hand, Drell-Yan production originates from valence and sea 
quark annihilation and hence the difference between the predictions of models A and B is small. 



fcoh 
'A 



dt 




t=0 



(56) 
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pA -» Xh + Y A 
g(pA)/Aa(pp) 



0.06 



pA -> Xh + A 
g(pA)/A 2 g(pp) (t=0) 




1 .channel 
model A 



0.1 0.2 0.3 0.4 0.5 



0.1 0.2 0.3 0.4 0.5 



Figure 8: The extra nuclear suppression factors Sa of (52) and (56) for the incoherent (pA — > 
XJi+Ya) and coherent (pA — > Xh+A) diffractive production of a system h carrying momentum 
fraction x in proton scattering on A = 64 (Cu) nuclei at corresponding proton-nucleon cm. 
collision energy of y/s = 300 GeV. The dashed, dotted and continuous curves correspond to 
using a single-channel and double-channel (models A and B) eikonals respectively. 



In Fig. 9 we include the prediction for the nuclear suppression in diffractive bb, as well as 
cc, production. The only difference is that the scale is larger: p, — 13 GeV rather than 5 GeV. 
By comparing cc and bb production we see the results are not sensitive to the choice of scale. 
In Fig. 10 we present the A dependence for the more realistic model B, for a large and a small 
value of x. 



4.3 Comments on proton— nuclear diffractive processes 

The above examples demonstrate the role that diffractive eigenstates (with different absorptive 
cross sections) play in the description of diffractive processes involving nuclei. In comparison 
with the over-simple single-channel Glauber approach, the nuclear suppression factors reveal 
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Figure 9: As in Fig. 8 but for proton collisions on A = 197 (Au) nuclei. 



an informative rich structure. They vary by more than a factor of two for the various hard 
diffractive processes. 

Figures 8-10 are qualitatively valid for all energies. The suppression factors of nuclear 
origin, Sa, decrease only gradually with increasing energy due to the slow growth of <j™ ot (s). 
Moreover, these nuclear shadowing effects for the interactions of protons with nuclei are similar 7 
to the suppression factors in proton-proton interactions at much higher energies [50]. (Recall 
that for both pA and pp collisions we refer to the proton-proton energy \/s.) For instance, we 
see from Fig. 8 that the nuclear suppression factor in p-Copper collisions is Sa ~ 0.1 at RHIC 
energies of yfs = 300 GeV, which is comparable to the suppression S p in hard diffractive pp 
collisions at the higher Tevatron energy, see for example, Fig. 6. This interesting correlation 
may be anticipated from the dependence of Sa of (51) and (52) on the exponential of the 

7 The similarity is not entirely complete as there is a difference in the impact parameter distributions for pA 
and pp collisions, and in the increase of the radius of the interaction with energy. 
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pA -» Xh + Y A P A -> Xh + A 




A 10 2 A 10 2 



Figure 10: The A dependence of the extra nuclear suppression factors Sa of (52) and (56) for 
the diffractive production of a system h in pA collisions. 

product — a™ ot (s)T A (b) . Total cross sections increase with energy as s A with A ~ 0.1, while 
Ta(P) ~ A 1 / 3 . Thus an increase of s by two orders of magnitude is equivalent to an increase in 
A by a factor ~ 4. Thus, since Sa for pA collisions with A ~ 64 at RHIC energies is comparable 
to S p for pp collisions at the Tevatron (s ~ (3-4) x 10 6 GeV 2 ), we expect that S p at the LHC 
(s ~ 2 x 10 8 GeV 2 ) will be close to S A for pA collisions with A ~ 197 (see Fig. 9). So the 
observation of shadowing of hard diffractive processes on heavy nuclei at RHIC may serve as a 
guide to the size of the suppression factors S p which occur in similar processes in pp collisions 
at the LHC energy. Of course, the S p factors can be calculated directly, but confirmation from 
a study of proton-nuclear scattering data at RHIC would be valuable. The S p factors need to 
be known, for example, in plans to identify New Physics phenomena in observations of such 
diffractive processes at the LHC. 
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5 Concluding remarks 



Investigation of both soft and hard diffractive processes provides important information on the 
dynamics of strong interactions at high energies. With increasing energy the role of unitarity 
effects (which, in the Reggeon approach, materialise as multi-Pomeron exchange interactions) 
becomes more and more important. The association of eigenstates of the diffractive part of 
the T-matrix with definite partonic configurations leads to extra insight into the dynamics 
of diffractive processes in QCD. It allows one to determine the effective sizes of the partonic 
configurations in hard diffractive processes. The use of nuclear targets provides new possibilities 
for detailed studies of various aspects of both soft and hard diffractive processes. We have 
demonstrated that interesting information on the properties of partonic configurations can 
be obtained from a study of hard diffractive processes in pA-collisions at RHIC and LHC. A 
knowledge of diffraction dynamics, and especially of shadowing effects in pp and pA collisions, is 
important for predictions of the expected manifestations of New Physics in diffractive processes 
at proton colliders, see, for example, [3] . A well known example is double-diffractive Higgs boson 
production, see, for example, [22, 4, 51]. 
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